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Abstract 

We establish a mortar boundary element scheme for hypersingular boundary integral equations 
representing elliptic boundary value problems in three dimensions. We prove almost quasi-optimal 
convergence of the scheme in broken Sobolev norms of order 1/2. Sub-domain decompositions can 
be geometrically non-conforming and meshes must be quasi- uniform only on sub-domains. Numerical 
results confirm the theory. 
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1 Introduction and model problem 

In the finite element framework, mortar methods are used to discretize a given problem independently 
on sub-domains. It is a non-overlapping domain decomposition method. Necessary continuity require- 
ments on interfaces of the sub-domains are implemented via Lagrangian multipliers. The motivation is 
to facilitate the construction of finite element meshes on complicated domains and to allow for paral- 
lelization. Bernardi, Maday and Patera introduced this technique and gave first analyses in ^[^. Later, 
geometrically non-conforming sub-domain decompositions and problems in R'^ have been studied by Ben 
Belgacem and Maday OH]- There is a large number of publications on mortar methods, all dealing with 
the discretization of differential equations of different types and with related numerical linear algebra. 
The first papers, just mentioned, derive a priori error estimates in the framework of non-conforming 
methods involving a Strang type estimate. 

In this paper we establish a mortar setting for the boundary element method (BEM) and prove almost 
quasi-optimal convergence for a model problem involving the hypersingular operator of the Laplacian. 
The advantages of this domain decomposition scheme (easier construction of meshes and availability of 
parallel techniques) also apply to the BEM. To be precise, we apply the mortar technique directly to 
the boundary element discretization, not as a coupling procedure between boundary and finite elements 
as in [9]. The analysis of finite elements for the discretization of boundary integral equations of the 
first kind goes back to Nedelec and Planchard [2T], and Hsiao and Wendland [17j. Stephan }22j studied 
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boundary elements for singular problems on open surfaces. Hypersingular boundary integral equations 
are well posed in fractional Sobolev spaces of order 1/2 and conforming Galerkin discretizations require 
continuous basis functions. Due to the non-existence of a trace operator in these Sobolev spaces, needed 
for the analysis of interface conditions, mortar boundary elements give rise to a variational crime. Indeed, 
it turns out that there is no well-defined continuous variational formulation of the mortar setting for the 
BEM. Instead we will analyze the discrete mortar scheme as a non-conforming method for the original 
un-decomposed integral equation. We follow the analysis presented in [1] where projection and extension 
operators are used to bound the approximation error in the kernel space (of functions satisfying the 
Lagrangian multiplier condition). Note that there is a shorter presentation by Braess, Dahmen and 
Wiener [7] where the simpler argument ^ Remark III. 4. 6] is used to bound this error by a standard 
approximation error (in un-restricted spaces). Nevertheless, in our case the Strang type error estimate 
has a more complicated structure and it is not straightforward to follow the argument |6i Remark III. 4. 6]. 

We will make use of some preliminary results in [131 I16j . In [13] we studied the discretization of 
hypersingular operators on open surfaces using functions that vanish only in a discrete weak sense on the 
boundary of the surface. Such functions in general do not belong to the energy space of the operator and 
require a different variational setting. This setting will be used also for the mortar boundary elements. 
In [16] this setting served to establish (non-conforming) Crouzeix-Raviart boundary elements and to 
prove their quasi-optimal convergence. Main tool in that paper is a discrete fractional-order Poincare- 
Friedrichs inequality. It serves to show ellipticity of the principal bilinear form of the discrete scheme. 
In this paper we generalize this inequality to the geometrically non-conforming case, needed for general 
mortar decompositions. Again, it is needed to prove (quasi-) ellipticity of the principal bilinear form. 
Our model problem is defined on an open flat surface T with polygonal boundary. We prove that, 
up to logarithmical terms, the mortar boundary element method converges quasi-optimally, subject to 
a compatibility condition of the boundary meshes and the meshes on the interfaces for the Lagrangian 
multipliers. Here we rely on the known Sobolev regularity of the exact solution leading to almost 0(/i^/^)- 
convergence where h is the maximum mesh size. Our techniques are applicable also to polyhedral surfaces 
and include meshes of shape-regular triangles and quadrilaterals. 

An overview of this paper is as follows. In the rest of this section we recall definitions of fractional 
order Sobolev norms and formulate the model problem. In Section [2] we define the mortar scheme 
and present the main result (Theorem 12. ip establishing almost quasi-optimal convergence of the mortar 
boundary element method. Technical details and proofs are given in Section [3l In Section H] we present 
some numerical results that underline the stated convergence of the mortar BEM. 

First let us briefly define the needed Sobolev spaces. We consider standard Sobolev spaces where the 
following norms are used: For a bounded domain S C R" and < s < 1 we define 



2 1 1 1 1 2 I 1 2 



with semi-norm 




For < s < 1 the space H'^{S) is defined as the completion of C^{S) under the norm 
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For s € (0,1/2), || • ||j5's(_5) and || • ||_f/=(S') are equivalent norms whereas for s € (1/2,1) there holds 

H^(S) = Hq(S), the latter space being the completion of C^{S) with norm in H'^(S). Also we note that 
functions from H^{S) are continuously extendable by zero onto a larger domain. For details see, e.g., 
[181 [Uj. For s > 0, the spaces H~'^{S) and H~'^{S) are the dual spaces of i/*(S') and H'^{S), respectively. 

Let r be a plane open surface with polygonal boundary. In the following we will identify T with a 
domain in R^, thus referring to sub-domains of F rather than sub-surfaces. The boundary of F is denoted 
by dr. 

Our model problem is: For a given f E -Z^^(F) find u E H^^'^{T) such that 



1 d f d 



Wuix) :=- — —/ uiy) — j—j dSy = /(x), x E F. (1.1) 



Here, n is a normal unit vector on F, e.g. n = (0,0,1)^. We note that W maps ^-"^/^(F) continuously 
onto i7^^/^(F) (see [23j). We have the following weak formulation of (jl.ip . Find u E i?^/^(F) such that 

{Wu, v)r = if, v)r E H'/\r). (1.2) 

Here, (•, •)r denotes the duality pairing between H-'^/^{T) and ^^/^(F). Throughout, this generic notation 
will be used for other dualities as well, the domain mentioned by the index. 

A standard boundary element method (BEM) for the approximate solution of (jl.2p is to select a 
piecewise polynomial subspace Xh C H^^'^iT) and to define an approximant Uh E Xh by 

{Wuh,v)r = {f,v)r e X^. 

Such a scheme is known to converge quasi-optimally in the energy norm, cf. Remark 12.11 below. In the 
numerical section we will compare such a conforming approximation with a mortar approximation, for 
the case where the meshes are globally conforming. 



2 Mortar method and main result 

In this section we introduce the mortar boundary element method for the approximate solution of the 
model problem (|1.2p . First we discuss the decomposition of F into sub-domains. Then we introduce the 
discrete approximation spaces. The main result of this paper is given at the end of this section. 



2.1 Sub-domain decomposition 

We consider a decomposition of F into non-intersecting sub-domains Fj, i = 1, . . . ,N, giving rise to a 
coarse mesh 

r := {Fi,...,Fjv}. 

For ease of presentation we assume that each Fj is either a triangle or quadrilateral. More general 
decompositions into polygonal sub-domains can be dealt with by further decomposing into triangles and 
quadrilaterals and by considering conforming interface conditions on additional interfaces. The mesh T 
can be non-conforming but must satisfy the assumption (Al) below. The diameter of a sub-domain Fj 
is denoted by Hi, and H := maxj=i^...^7v -f^i- The interface between two neighboring sub-domains Fj, Tj 
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{i j, Fj nFj contains more than a point) is denoted by 7ij. For our analysis below we need the following 
assumption. 

(Al) Each interface 'jij consists of an entire edge of Fj or Tj. 

The (relatively) open edges of a sub-domain Fj are 7^, j = 1, . . . ,m. Here, m is a generic number 
[m = 3 if Fj is a triangle and m = 4 otherwise). Using the symbol dT for the boundary of F, and similarly 
dTi for the boundary of Fj, the skeleton of the sub-domain decomposition is 

7:= ul,dTi\dT. 

According to assumption (Al) the skeleton is covered by a set of non-intersecting edges ^ij. We number 
the edges like 71, . . . , 72,, giving a decomposition of the skeleton like 7 = U^^t";. In the following we will 
denote this decomposition of the skeleton by 

r := {71,..., 7l}. 

We will refer to these edges as interface edges. Each interface edge 7; is the interface between two sub- 
domains Fj, Fj and is an entire edge of one or both of them. Given an integer / (1 < / < L) we denote 
by ^lag (respectively, /mor) the number of a sub-domain which has 7; as an edge (respectively, the number 
of the other sub-domain). 

As mentioned before, the selection of the index pair (Ziag, ^mor) for / € {1, ... , L} is not unique but will 
be fixed for a specific sub-domain decomposition of F. Below, we will introduce a Lagrangian multiplier 
on the interfaces and on ji we will use a mesh related to the mesh on F^j^^ . The side of 7/ stemming from 
^imor is usually called mortar side in the finite element literature and this explains our notation. The side 
defining the Lagrangian multiplier is often called non-mortar side. 

Corresponding to the decomposition of F we will need the product Sobolev space 

N 

H%T) := II H'{K) = l[H'{Ti) 
K&r i=i 

with usual product norm. 

2.2 Meshes and discrete spaces 

On each of the sub-domains Fj (i G {1, . . . , N}) we consider a (sequence of) regular, quasi-uniform meshes 
7j consisting of shape regular triangles or quadrilaterals, Fj = U^gT-.T. The maximum diameter of the 
elements of % is denoted by hi and we use the symbols 

h := min hi, h := max hi. 

~ i=l,...,N 1=1,. ..,N 

Throughout the paper we assume without loss of generality that h < 1. This makes the writing of 
logarithmic terms in h easier. 

In the case of F being a square. Figure 12.11 shows a conforming sub-domain decomposition (a) and a 
non-conforming sub-domain decomposition (b), both with globally non-conforming meshes. 
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Figure 2.1: Sub-domain decompositions with non-conforming meshes. 

Now we introduce discrete spaces on sub-domains consisting of piecewise (bi)hnear functions, 
Xh,i '■= {v £ C^{yi)'t v\t is a polynomial of degree one VT G %, ulamari =0}) i = 1, . . . ,N. 
The global discrete space on F is 

TV 
1=1 

Note that functions v € do satisfy the homogeneous boundary condition along dT but are in general 
discontinuous across interfaces. Therefore, Xh is not a subspacc of the energy space H^/^(T). Functions 
from different sub-domains will be coupled via a discrete Lagrangian multiplier on the skeleton. To this 
end we introduce a mesh on the skeleton 7 as follows. 

On each interface edge 7/ there is a trace mesh 7/, L, inherited from the mesh 7/, on the sub-domain 

^ ' * f-lag I It ''lag 

F/j^^. (We recall that by definition, 7/ is an entire edge of the sub-domain with number /lag-) This trace 
mesh is quasi-uniform with mesh width hi^^^. Now we introduce a new (coarser) quasi-uniform mesh Qi 
on 7/ in such a way that the following assumption is satisfied. 

(A2) For any Z G {1, . . . , L} there holds: the mesh Qi is a strict coarsening of the trace mesh 'Tl^^JJ^. In 
particular, any interior node of Ti^^J^^ together with its two neighboring elements (intervals) is covered 
by one element of Qi. 

The mesh width (length of longest element) of Qi is denoted by ki, and k := max;=i fc;. On each 
interface edge we define a space of piecewise constant functions, 

Affcj := {v G L?{'^i); v\j is constant VJ € Qi], I = 1, . . . , L. 

The space for the discrete Lagrangian multiplier then is 

L 

Mk ■■=Y[Mk,i. 

1=1 



5 



Notations. The symbols "<" and ">" will be used in the usual sense. In short, ah{v) < bh{v) when 
there exists a constant C > independent of the discretization parameter h and the involved function 
V such that ah{v) < Cbh{v) for any v of the given set. The double inequality ah{v) < bh{v) < ah{v) 
is simplified to a^(u) ~ bh{v). The generic constant C above is usually also independent of appearing 
fractional Sobolev indexes e > 0, but this will be mentioned. We note that these notations usually do 
not mean independence of involved constants on the decomposition T of T. In this paper we consider a 
generic decomposition T which is fixed and estimates will in general depend on T. 

Throughout the paper we also will use the notation Vj for the restriction of a function v to the 
sub-domain Tj. 

2.3 Setting of the mortar boundary element method and main result 

For the setup of the mortar boundary element method we need some operators. We introduce the surface 
differential operators 

curl if := {dx2^P, -dx^f, O) , curl := dx^f2 - dx^^i for (f = {ipi,ip2, ^fs)- 

The definitions of the surface curl operators are appropriate just for flat surfaces (as in our case) but 
can be extended to open and closed Lipschitz surfaces, cf. [Sj |T3]. We define corresponding piecewise 
differential operators curl// v and curl// (p by 

N N 

curl// u := ^ (curlr, Vif , curl// </> := ^ (curlr, . 

i=l i=l 

The notations curlr, and curlr^ refer to the restrictions of curl and curl, respectively, onto Ti, and {■)^ 
indicates extension by zero to F. We made use of the notation introduced before, Vi = f Ir^, (pi = (plrt- 
Furthermore, we need the single layer potential operator V defined by 

V<p{x) [ dSy, <p € (^-i/^(F))3, X G F. 

Air Jy \x-y\ 

For the formulation of the mortar boundary element method we define, for sufficiently smooth functions 
V, w, fj,, the bilinear forms a(-, •) and •) by 

N 

a{v,w) := (y curl// curl// t(;)7- := ^^(y curl// t), curlpj tL')ri, 

1=1 

L 

b{v,fi) := {[v],n)r ■■=^{[v],n)-yi. 

1=1 

Here, as mentioned before, for a domain 5 C F or an arc S, (•, ■)s denotes the L^(S')-inner product and 
its extension by duality, and [v] is the jump of v across 7, more precisely 

M Ul = ^«lag \ll - Vlmor I7P l = l,...,L. 

Of course, for sufficiently smooth functions v, w, fj, there holds 

a{v,w) = curl// curl// tt;)r, b{v,^) = {[v],fj,)^. 
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Note that we will use the introduced notations (•, •)r and (•, •)t for duality pairings of product spaces 
corresponding to the given decompositions (T and r). We also define, for a sufficiently smooth function 
V, the linear form 

N 
i=l 

where / G -^^^(r) is the function given in (jl.ip . The mortar boundary element method for the approximate 
solution of (jl.2p then reads: Find € and € such that 

a{uh,v) +b{v,\k) = F{v) \/v £ Xh, n 

biuh,i^) =0 VV^GMfc. ^ ■ ^ 

This scheme is equivalent to: Find such that 

a{uh,v) = F{v) yv € Vh 

where 

Vh = {ve Xh-, b{v, ij)=0 VV G Mk}. (2.2) 
The main result of this paper is as follows. 

Theorem 2.1. There exists a unique solution {uh,Xk) of ()2.ip . Assume that the solution u of ()1.2p 
satisfies u G ^^/2+^(r) (r € (0, 1/2]). Then there holds 

\\u-Uh\\H^/2(^r) ^ dlog/iP^'' + |l0gZlP''^^'^)ll^^lli/l/2+r(r). 
For proportional mesh sizes h and k this means that 

\W - Uh\\Hi/2(j^ < I log/^p/l'■||u||^l/2+r■(^)• 
T/le appearing constants in the estimates above are independent of h and k provided that the assumptions 
on the meshes, in particular (Al) and (A2), are satisfied. 

A proof of this theorem is given at the end of Section [3l 

Remark 2.1. As in [13j we note that, in our case of an open surface T, the solution u of (jl.ip has strong 
corner and corner-edge singularities which cannot be exactly described by standard Sobolev regularity. It 
is well known that u G H^{T) for any s < 1 (see, e.g., [24JJ so that the error estimate by Theorem 12.11 
holds for any r < 1/2. In general u Hq{T) but a more specific error analysis for the conforming BEM 
yields for quasi-uniform meshes the optimal error estimate 

\W-Uh\\Hi/2(j^) < /^^''^ 

see [5]. The logarithmical perturbations in h of our error estimate are due to the non- conformity of the 
mortar method. They stem from the non-existence of a trace operator within II^^'^(r) and from non-local 
properties of the fractional order Sobolev norms (the difference between H^/'^ and H^/'^ -spaces). 
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3 Technical details and proof of the main result 



We start by citing some technical results (Lemmas 13. m3.4|) which are needed to deal with the fractional 
order Sobolev norms. Afterwards we study a discrete Poincare-Friedrichs inequality (Proposition 13. ip 
which will be applied to prove ellipticity of the bilinear form a(-, •) on Vh- Afterwards an integration-by- 
parts formula for the hypersingular operator is recalled from |13j and adapted to our situation of many 
sub-domains. Then, Lemma [331 states the well-posedness of integration by parts. Lemmas 13.6113. 131 study 
requirements for the Babuska-Brezzi theory and provide details for a Strang-type error estimate which is 
given by Theorem 13. 11 Later, Lemmas 13. 141 and 13. 151 are needed to analyze the bound of the Strang-type 
estimate and lead to Theorem 13.21 which gives a general a priori error estimate for the mortar BEM. The 
section is finished by giving a proof of the main result (Theorem 12. ip . 

Lemma 3.1. [15] Lemma 5] Let R C be a Lipschitz domain. There exists C > such that 

Mhhr) < 1/2- VsG (-1/2,1/2), yveH^R). 

Lemma 3.2. Let R C be a Lipschitz domain, and let v be a piecewise linear function defined on a 
quasi-uniform mesh on R with mesh size h < 1. There exists a constant C > which is independent of 
h (but may depend on R) such that there holds 

\Hh-1/2^r) < C\logh\ \\v\\jj-i/2^jiy 

Proof. By Lemma 6] there holds for a piecewise polynomial function of degree p the estimate 

p + 1 

||t'||j^-i/2(j^) < Clog{—^)\\v\\jj-i/2^ji-j, p>0,h<l. 

Fixing p gives the claimed bound. The proof of [15, Lemma 6] gives full details for rectangular meshes. 
For triangular meshes the proof applies as well by making use of Schmidt's inequality for triangles, cf. [IH 
Lemma 5.1]. Nevertheless, we are considering only polynomials of low degrees where Schmidt's inequality 
is not needed. □ 

Lemma 3.3. |13l Lemma 4.3] Let i? C R^ be a bounded Lipschitz domain. There exists C > such that, 
for any e G (0, 1/2), there holds 

(J 

\HLHdR) < ^\Mm/2+.(R) yv G H^I^+'{R). 
Here dR is the boundary of R. 

Lemma 3.4. For S being one of the sub-domains Pj G T or F there holds 

\v\hi/2(^S) ^ il curl5 v\\h;^/^^s) ^ H'^\S). (3.1) 
The restriction of curls onto H^/'^{S) is continuous, 

curls: H^I\S) ^ H;^'\s). (3.2) 

Moreover, there holds the continuity 

curls : ^^/^+"(S) ^ H-'^'^^\S) Vs G [0, 1/2]. (3.3) 
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Proof. The bounds (|3.ip and (j3.2p are proved by Lemmas 4.1 and 2.2 in [13], respectively. By Lemma 2.1 
in [13], curls : H^/^{S) H~^''^{S) is continuous, and curl^ : H^{S) L^iS) = H^{S) is continuous 
as well. Estimate ()3.3p then follows by interpolation. □ 



The following result is a generalized version of a discrete Poincare-Friedrichs inequality in fractional 
order Sobolev spaces, cf. Theorem 8 in [16j . 

Proposition 3.1. There exists a constant C > 0, independent of the decomposition T as long as sub- 
domains are shape-regular, such that for all e € (0,1/2] there holds 

Here, \^i\ denotes the length ofji. 

Proof. For the case of conforming decompositions T of F into triangles, [16^ Theorem 8] proves that there 
holds 



Hh^r) < C (^e">lW(r) + El^'l"'"''(/ Md.)' + \l^vdx\'^ yv G H'/'+^{f). (3.4) 



It is easy to see that the mean zero term can be avoided by assuming the homogeneous boundary 
condition for v. To obtain the result for our non-conforming decomposition T including quadrilaterals 
we introduce further edges to reduce quadrilateral sub-domains to triangles and to transform T into a 
conforming decomposition T. By definition of the Sobolev-Slobodeckij semi-norm there holds 



-/2+^(r) - l^!/fi/2+^(r)- 



We note that for new edges 7' there holds [v] \y = by the trace theorem and the regularity v G H^/'^~^''{T) 
for any T ^T. The result then follows from (j3.4p . □ 



Following [13] we now examine an integration-by-parts formula for the hypersingular operator. For a 
smooth scalar function v and a smooth tangential vector field Lp, integration by parts gives 

(curlr, V, (^)r, = {v, curlr, (^)r, - (u, • ti)ar, , i = 1, . . . , iV. 

Here, tj is the unit tangential vector on dVi (oriented mathematically positive when identifying Fj with 
a subset of which is compatible with the identification of F as a subset of R^). Applying this formula 
to (/? = (y curlru)[ri; we obtain for smooth functions v and u 

{ti ■ V curlr n, Vi)dT^ = (curlr, V curlr n, Vi)r, - {V curlr n, curlr, Vi)r, ■ 

Now we sum over i and take into account that ti = —tj on jij. Further we let 70 := dT, use the convention 
for the jump [v]\^^^ = v\.yg, denote by to the unit tangential vector along dT (again mathematically positive 
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oriented) and let Oiag := (remember the notation /lag and /mor for the numbers of the Lagrangian 
multiplier side and mortar side of 7/, respectively). This yields 

L N N 

^(t/j^g • y curlr n, [w])^, = (curlr, V curlr u, Vi)ri - y^(y curlr tt, curlr, Vi)ri 

1=0 i=l i=l 

= (curlr V curlr u, v)j- — {V curlr u, curl// v)j- 

for a piecewise (with respect to T) smooth function f on F with Vi := vIy^, as defined before. In the last 
step we used the fact that 

curl// til = curlr If \/w € Hy^{r), 

which holds by a density argument and the continuity of curlr : Hl^'^iT) H~^/'^{T) as the adjoint 

operator of curlr : H^^^iT) H^^^^{r), cf. (fO) . 
Now we use the relation 

Wu = curlr V curlr u (n G H^^^{T)), 
see [191 120] and [131 Lemma 2.3]. Then choosing a piecewise smooth function v with tijgr = we obtain 

L 

(A,H), = ^(A,H)^, = {Wu,v)t- {V curlr u,curlHv)T. (3.5) 

1=1 

Here, A denotes our Lagrangian multiplier on the skeleton 7 defined by 

Myi ■=h,g- curlr u)\^^^, l = l,...,L. (3.6) 

Relation (j3.5p does not extend to u G H^^'^{T) since the trace of such a function v onto 7 is not well 
defined. However, there holds the following lemma. 

Lemma 3.5. Forn G H^I'^{T) with Wu = f e L^{T), ([33]) defines A G ULi H^'ili) for any s G (0, 1/2]. 

Remark 3.1. The above lemma can be extended to values of s larger than 1/2. Though small values of 
s represent the interesting cases, the limit s = being excluded. Also, the condition on f can be relaxed 
but excluding the case f G H^^/^(r) which is the standard regularity using the mapping properties of the 
hypersingular operator. 

Proof of Lemma 13.51 We must show that A defined by (j3.5p is a bounded linear functional on 
Uti H'ili), the dual space of Uti ^"'(tO- 

Let V G Y[h=i -^'^(li) t'G given. We continuously extend v to an element v G H^^^/'^{T) with v = Q 
on dV such that [i;]|^, = vj^^. (Simply extend v on each interface edge 7/ to a function in H^'^^/'^iVi^^ ) 
vanishing on dVi^^ \ji and extend by zero to the rest of T. Then sum up with respect to L) The definition 
of A is independent of the particular extension v, see |13j for details in the case of one sub-domain. Using 
a duality estimate we obtain from (13. Sp 

L 

X]('^iH)7i = (/,S)r - curlr n, curl// ■u)r 
1=1 

N 

< ll/llL2(r)ll^^llL2(r) + ^ |iycurlru||^i/2--(P^)|| curlr, Vi||^.-i/2(P^^. (3.7) 
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Now, for s G (0,1/2] the norms in -ff^^ '^{^i) aiid -ff^^ '^{^i) are equivalent (cf., e.g., [15]) so that 
together with the mapping property of V [10] , 



and (j3.2p we obtain 

Af AT 

E||ycurlru|P. i/2_. < ||y curlr ulp i/2_s. . 5. II curlr tilp i/o-s, . 5, ll'uli m /9^T-\- 
II i (Fi) ^ "Mt'^ (Fi) ~ " ^ "Mt-^^ ^(P) II ll_H-i/2(r) 

i=l i=l 

Here, the appearing constants are independent of u but may depend on s. Also, using (j3.3p we are able 
to bound (with constant independent of v) 



2 



N N 



J2 II curlr, i)i||^,_i/2 < J] ||vi||^,+i/,,^ ^ = \\v 



i=l i=l 

Taking the last two estimates into account, p.7p proves that 

L 



X]^-^'M)7i ~ (ll/llL2(r) + ||n||jji/2(r)) WvUh^+^/^t)- 



1=1 

Using the continuity of the extension (with constant independent of v) 

L 

||l'||^s+l/2(^) < ^ ll^ll^=(7i) 
1=1 

finishes the proof. □ 
Lemma 3.6. 

||curl^^u||2^_i/2^^^ > Ilog/ir^|l?;|t^i/2(^) Vf £ Vh 

Proof. By (jS.ip there holds 

N N 

||curl/^w||^_i/2^^^ > ll^"''lr,^;i||^-i/2(^^^ > \^^\m/Hr,) = I^Ihi/2(t) ^ H^^\T). (3.8) 

i=l j=l 

For V € Vh there holds / [u] ds = for any interface edge 7; since by construction contains the 
piecewise constant function which has the value 1 on 7; and vanishes on 7 \ 7^. For the definition of Vh 
(the discrete kernel of b{-, •)) see (|2.2p . Therefore, Proposition 13.11 proves that 

ll^llL2(r) ~ e"^l'^l^i/2+.(r) e ^h- 
Here, the appearing constant is independent of e € (0, 1/2]. Making use of the inverse property we bound 

N N 



i=l i=l 
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so that, with the previous estimate, 

Iblli2(r) < e"'r'l^^lW(r) V^; e y,.. (3.9) 

Selecting e = | logh\~^ (for h being smaU enough) and combining ()3.9p with (13. 8p proves the statement. 

□ 

Lemma 3.7. The bilinear form a(-, •) is almost uniformly Vh-elliptic. More precisely there hold the lower 
hounds 

a{v,v) > |log/ir^||v||^i/2(^) Vt; € Vh 

and 

a{v,v) >\logh\'^/'^\\v\\jji/2^r)\\curlHv\\^-i/2^^^ G 14. (3.10) 
Proof. First we note that for v ^ Vh C L'^(r) there holds ycurl/^-f G L'^(r) so that 

{V curljy u, curl// f )7- = {V curl// w, curl// v)r- 

Using the ellipticity of V : H^^ ^ (F) — > /^^^(F) and Lemma 13.61 we then obtain for v Vh 
a{v,v) = (1/ curl// curl// i;)r > || curl// _i/2^^^ > | log ||u||^i/2(r)' 

which is the first assertion. The estimate (|3.10p is obtained by bounding || curl// 1/2^^^ only once 
with the help of Lemma 13.61 □ 

Lemma 3.8. The bilinear form a(-, •) is almost uniformly continuous on Xh- More precisely there holds 

a{v,w) < \ loghf \\v\\jji/2^r)\\M\m/2{r) ^^'^ ^ ^h- 

Proof. As in the proof of Lemma 13.71 we note that for v, w ^ Xh C ^^(F) there holds 

{V curl// V, curl// w)t = {V curl// v, curl// w)r. 

— 1/2 1/2 

Then, using the continuity of V : iJ^ (F) — > H^^ (F) and the estimate for fractional order Sobolev 
norms || • || --i/2_, < || • || --i/2,_,, we obtain for v,w e Xh 

a{v, w) = {V curl// v, curl// w)r < || curl// v\\ ^,-1/2.^. \\ curl// f^-1/2.^. 

< ||curl//t;||^_i/2^^^||curl//u;||^_i/2^^^. (3.11) 

Now making use of Lemma 13.21 and p.3p we bound 

llcurlr, Uif i/2_^ < |log/iip||curlr, fif 1/2,^, < | log /iip||ui||^i/2fr )' 

giving 

II curl// ?;|[ ^^-i/2(^^ < I log /i| ||f ||//i/2(r) G X/^. 
Combination with (j3.1ip proves the statement. □ 
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For the eventual error estimate we need the boundedness of the bilinear form a(-,-). However, 
Lemma 13.81 is not applicable to non-discrete functions. Instead we will use the next lemma. 

Lemma 3.9. Assume that u G H^/'^^'^iT) (r > Oj. Then there holds 

a{u -v,w) < s^^\\u - ■u||^i/2+s(7-)|| curl/^ ?i;||^-i/2^^^ yv,w G Xh, Vs G (0,min{r, 1/2}]. 

In particular, the appearing constant is independent of s. 

Proof. First we note that (j3.1ip holds also for continuous functions, so that 

a{u -v,w) < ||curlH(^i - || - -1/2 || curlj/ - -1/2 yv,w G Xh. 

Using the continuous injection H^^ ^ (Fj) ^ (Fj) and Lemma |3. II we bound for i G {1, . . . , N} 

\\c\yr\Y^{ui-Vi)\\^-i/2^^^^ < \\c\iT\Y^{ui-Vi)\\^-i,2+s^^^^ < s'^^Wcurlr^iui - Vi)\\^-i/2+s^^^y 

The continuity of curlr, : H^/^+'{Ti) H~^^'^~^' (Ti) for any i G {1,...,A^} by finishes the 

proof. □ 

In order to analyze the error bound of the Strang-type estimate by Theorem 13.11 below, we need to 
extend functions from interface edges to sub-domains. This is also required to prove an inf-sup condition 
for the bilinear form b{-,-). 

To this end let us define extension operators that extend piecewise linear functions from interface 
edges to piecewise (bi)linear functions on the corresponding Lagrangian sub-domain. 

El : Xh,i^^J^i ^/i,/iag' l = l,---,L. (3.12) 

Here, for v G X^^ii^ I7, the extension Eiv is defined as the function of X^^i^^ that coincides with v in the 
nodes on 7; stemming from the mesh T/j^^ and is zero in the remaining nodes of Ti^^^ . 

Lemma 3.10. 

\\EM\H^(Ti^^^)<h]l^~"\\v\\L2(^^^^ Vt; G Xfc,,,^J^,, VsG[0,l], Z = 1,...,L. 
In particular, the appearing constant is independent of s. 

Proof. Using the equivalence of norms in finite dimensional spaces and scaling properties of the L^-norm 
one obtains, by taking into account the construction of Ei, 

\\EM\h{T,^^^) ^ KJv\\l2(^^^) Vt- G Xh,i,J^,. 

Analogously we find 

\\Eiv\\]mr, ) = ll-^/^lli2(r ) + |Si7;|^i(p ) 



lag l7i • 



The result then follows by interpolation. □ 
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Lemma 3.11. The bilinear form satisfies the discrete inf-sup condition 



3/3 > : sup ,, > /? ||^||l2(^) V/x G M^. 

^e^A{o} IPll//i/2(r) 

Here, the constant (3 is independent of h and k subject to the assumptions made on the meshes. 

Proof. Let /i € be given. On each interface edge 7^, /x is a piecewise constant function on a mesh 
that is coarser than the trace mesh Ti-^^J\^i stemming from the Lagrangian side L/j^^, cf. assumption (A2). 
On 7; we construct a piecewise hnear function wi G ^h,ii^J\'yi in the following way. For each element 
J ^ Qi, wi vanishes at the endpoints of J, coincides with /x at one interior node of J and is linearly 
interpolated elsewhere on 7; . See Figure 13.11 for an example where ^ is represented by the dashed line 
and wi by the solid line. The bullets indicate the nodes of the mesh for the Lagrangian multiplier and 
the dashes indicate additional nodes of the trace mesh (from the Lagrangian multiplier side). 

We then extend wi to wi in by first extending to Eiwi G X/^^/j^^, cf. (13.12p . and then further by 

zero onto F. Eventually we define v := X^^x 

Note that wi vanishes on all interface edges except 7;. The trace of wi onto 7; from F^j^ equals wi 
whereas the trace coming from the other side Fj^^^^j. vanishes. This yields 

[v] = [wi] = wi on 7i, I = l,...,L. (3.13) 

By the construction of wi there holds, uniformly for /i G M^, 

M\h(^ii)~ {wu^j)^i~\\wi\\\2(^^^^y l = l,...,L. (3.14) 

Also, taking into account that each sub-domain Fj has a limited number of (interface) edges, determined 
by the relation I £ {1, . . . , L} : /lag = i, Lemma [3.101 yields 



l7;l|2 

l^lli/i/2(r) 



^ 2 



E II E 

i=l ie{l,...,L}: iiag=i 



N L L 

< E ii^'«^HiW(r,) = Eii^'^'iiW, )<EiNiii^(7,) (3-1^) 

i=l iG{l,...,L}: iiag=i '=1 '''^ '=1 



Now, using (|3.13p . (|3.14p and ()3.15p . we finish the proof by bounding 
L L L 



1 /2 

Kv,fJ') =Y{[v],fJ')'yi =Y{'Wl,fJ^)-yi ^ II/^IIl2(7)(E 11^' 11^2(7;)) ~ \\f^\\L2{'y)\\v\\m/^T)- 

1=1 1=1 1=1 

□ 

Lemma 3.12. 

IIMIli2(7) ^ \ ^ogh\ \\v\\1j,/2(^r) 



14 



Figure 3.1: Construction of wi in the proof of Lemma 13.111 

Proof. By the triangle inequahty and Lemma 13.31 there holds uniformly for e G (0, 1/2) 
L 

IIHIli2(^) < '^{\\Vll^J\h{'yi) + \\'"lmor\\h{-yi)^ 
1 = 1 

- ^"^l^ ( ll^'lagll^l/2 + ^(r,. ) + ll^'morll^l/2+.(r,^^j) <e"i^^llHl/2 + .(r) ^^GXh- 
1 = 1 \ ''B / 

The inverse property, applied separately to Vi = v\r^, yields 

and selecting e = | log (for /i being small enough) finishes the proof. □ 
Lemma 3.13. The bilinear form b{-,-) is almost uniformly discretely continuous, in the sense that 

^(■W)/") ^ I logZil^^^ lbl|//i/2(r)ll/^llL2(7) Vt! € X/i, V// G Mfc. (3.16) 
Moreover, for given tp G L^(7), there holds 

b{v,i^) <\'^ogh\ inf IIV' - /x||i2( )|| curl/^-vll --1/2 Vu G T4. (3.17) 

Proof. There holds 

L 

Hv,i^) = Y,{[v],i^)^, < IIHIlL2(7)ll/i||L2(^) v^; G Xh,y^^ g m^. 

Estimate (j3.16p follows with the help of Lemma 13.121 To prove (j3.17p we start as before and note that 
by definition of Vh there holds b{v,ii) = V^u G M^, G Vh- Therefore, for any /i G and f G 14 we 
find that 

b{v,tl;) < ||H||i2(^)||V'-/i||L2(^). (3.18) 
The proof of (|3.17p is finished by noting that combination of Lemmas 13.121 and 13.61 yields 

IIHIIl2(7) < \logh\^/^\\v\\Hi/2^r} ^ \ ^ogh\ ||curlHw||^-i/2^^^ Vv G T4. 

□ 
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We are now ready to prove the following Strang-type error estimate. 



Theorem 3.1. System (j2.1|) is uniquely solvable. Let u and Uh be the solutions of (jl.2|) and (|2.ip . 
respectively. Assuming that u G H^/'^^'^iT) (r € (0,l/2]j there holds 



\u-Uh\\m/2(T)^\'^ogh\^^'^ is iid \\u-v\\Hi/2+s(r)+ sup 



uniformly for s £ (0, min{l/2, r}] . 

Proof. The existence and uniqueness of (u/j, Ajt) G x follows from the Babuska-Brezzi theory. 
Indeed, the bilinear form a(-,-) is continuous on Xh by Lemma 13.81 and V/i-elliptic by Lemma 13.71 and 
the bilinear form h(-, •) is continuous on X^ x by (j3.16p and satisfies a discrete inf-sup condition by 
Lemma 13.111 The continuity and ellipticity bounds depend on h but that does not influence the unique 
solvability of the discrete scheme. 

The error estimate is obtained by the usual steps. Combining the triangle inequality, the non-standard 
ellipticity and continuity properties of a (•,•), cf. ()3.10p and Lemma \3M we obtain for any v £ Vh 

\W - Uh\\m/2(T) ^ 11^ - f^ll/fi/2(r) + lb - ^/ill/fi/2(r) 



< 



I II 1 M f,|i/2 aiv-Uh,w) 

\u-v\\Hl/2rr) + \logh\ ' sup 7 

^ ' w€Vh\{0} CUtIhW .-1/2 



H, ' (r) 



1,11/2/ a{v-u,w) a{u-Uh,w) 
< If - ■y||/i'i/2(r) + I logiiil I sup II , — — h sup 



\w<^Vi,\{0} II curl// weVh\{0} II curl// w\\ ^-1/2^^^ J 

^ II II -In ,i1/2ii 11 11 ,i1/2 CL(u — Uh,w) 

< \\u-v\\fji/2(r) + s \logh\'\\u-v\fji/2+s(r) + \iogh\' sup 7 . 

^ ' ^ l^) «,ev,\{0} l|curl//ii;||^-i/2^^^ 

This proves the stated error bound. □ 

In order to analyze the upper bound provided by Theorem 13.11 we need, apart from the extension 
operators Ei defined before, projection operators tti acting on L'^i'yi) and mapping onto special continuous, 
piecewise linear functions on 7/, / = 1, . . . , L. We recall that on each 7; we have two meshes: the trace 
mesh '?ij^g|7, stemming from the mesh on the sub-domain F^j^^ of the Lagrangian side, and the mesh Qi 
for the Lagrangian multiplier. For each element J £ Qi we consider a hat function (pij that vanishes at 
the endpoints of J and has the tip at a node of Ti^^ that is interior to J. This choice is not unique if 
J contains more than two elements of the trace mesh. In that case we select an arbitrary but fixed node 
for the definition of (pi.j. Using this notation we define 

TTi : L'^i'yi) span{0;,j; J £ Gi} C ^Mi.gb,, l = l,...,L, (3.19) 

such that the integral mean zero conditions 

{v-TTiv,i)j = yj eGi, l = l,...,L, 

hold. This operator satisfies the following properties. 
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Lemma 3.14. For any v E L (ji), ttiv vanishes at the endpoints of ji, I = 1, . . . , L, and there holds 

{v-7riv,fi)^,=0 yveL\ji),yneMk^i, l = l,...,L, (3.20) 

hMlL^-ri] ^ II^IIl2(70 e L^M, l = l,...,L. (3.21) 

Proof. For / E {1, . . . , L} let u G L'^{'~fi) be given. By definition of tti, ttiv vanishes at the endpoints of 7/, 
and the orthogonality (j3.20p follows by noting that any E M^^i is constant on any J € Gi- 

To show (j3.2ip let J E ^/ be given. With cpij being the hat function defined previously (with height 
1) there holds 



TTlv = |^( j vds) (1)1 J on J 



so that 



Summing over J Gi finishes the proof. □ 

We are now ready to analyze the first term of the upper bound provided by Theorem 13.11 
Lemma 3.15. For r E (0, 1/2] let u E H^/'^+''{T). There holds 

L 

mf^ \\u - ^;|l^v2+.(r) ~ 11^ - M\li/2+s(^r) + Yl ^Ta^ ("^ " "^hjlh^ + 11^ " win^oMh^) ^ 
uniformly for s E (0, min{l/2, r}] . 



1=1 



Proof. Let w E Xh be given. We adapt w such that the new function satisfies the jump conditions 
defining Vh, cf. 1^. We set 

L 

V := w 

1=1 

with 

I elsewhere. 

Here, Ei and tti are the extension and projection operators specified in (j3.12p and (|3.19p . respectively. 
Note that, since T^iiwii^J-yi — ^/mork;) vanishes at the endpoints of 7/, the extension El^Tl{wl^^J■y^ — ui^^^J^J 
vanishes on dTi^^^ \ 7^. Therefore, using (j3.20p one obtains 

= (^ilag - + TT^^^/lag \ll - ItJ' /^)7i =0 ^ Mk,h l = l,...,L. 

That is, V £ Vh- We start bounding the error by 

N 



^1/2+^(7-) — Y2 |p* ^* ^ 

j=l /e{l,...,L}: /iag=i 



2 

//i/2+=(r,) 



AT L 

^ X] ll""* ~ '"^*llHi/2+^(ri) + X] ll'''ll/fi/2+«(r,j^^)- (3.22) 



i=l 1=1 
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Applying Lemma 13.101 (j3.2ip and the triangle inequality we find that there holds 

^ ^/Tag ( 11^ -^'lag 11^-^(70 + 11^ 11^^2(70) ■ (3-23) 

Combining (j3.22p and (j3.23p one obtains the assertion. □ 

The next result provides an a priori error estimate for the mortar BEM. 

Theorem 3.2. Let u and he the solutions of ()1.2p and ()2.ip . respectively. Assuming that u € 
(V S (0, l/2]j there holds A G 11^1-^^(70 ^''^^ have the a priori error estimate 

\\u-uh\\]j^/2^^r) ~ s-'^\^ogh\ i^\u-v\\]j^/2+,^j^ + h-^'Yl{\\^-'"hJ\^^^ 

+ I logZil^ ||A - /i||i2(^) Vf € Xh, V^i G Mfc 
uniformly for s E (0,r]. Here, A is t/ie Lagrangian multiplier defined by (j3.5p . (j3.6p . 
Proof. Since ti G i7^/^+''(r) there holds 

=t«iag • ("^curlru)!^, (^H'^M, l = l,...,L. 

To this end note that curlr : H^/^^^'{T) h[ ^^^(T) (combine (j3.2p with the continuity curlr : 
i7Q(r) — > ^^(r)) and V : ^^^(r) ^ ^1^'^^^ (T). The trace theorem concludes the claimed regularity 
of A. In particular there holds A G L'^i^)- 

By definition of V^, and making us of Lemma 13.5^ we find 

a(n — Uh, w) = a{u, w) — F{w) = —b{w, A) \/w E Vh. 

Application of (j3.17p yields 

a{u - Uh,w) <\logh\ inf IIA - ^||i2( )|| curl//?i'||£v-i/2,^, Vui e V^. 

Therefore, combining Theorem 13. II with Lemma 13.151 we obtain 

\\u - uh\\l^/2^r) ^ |log/i||s"' (^ll^-^llW+={r)+Z]^C(ll''"''''-JI''(70 + ll^-^'-^lli^(7o)j 

+ I log/ip ||A - ^Ili2(^)} V^; G X;,, V/x G M^. 

This proves the statement. □ 

Proof of Theorem 12. IL By Theorem 13.11 system (12. ip is uniquely solvable. We employ the general a 
priori estimate by Theorem 13.21 to show the given error bound. By standard approximation theory there 
exist V G Xh and n G such that 



k-^^ll^i/2+.(r) ^^^^'^ '^ll^ll|i/2+.(r) and ||A - /xHi^^^^ < A;^'^ ^ ||A||2^.(^^), 

1=1 
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and as in the proof of Theorem 13.21 one concludes that X^^i 11-^11/^(7,) ^ ll'"ll^i/2+r(p)- Lemma [331 
one bounds 

and accordingly the mortar part \\u — ■vZmor 11x2(7,)- Using these bounds in Theorem 13.21 and selecting 
s = I logh\^^ one obtains the assertion. □ 



4 Numerical results 

We consider the model problem (jl.2p with T = (0, 1) x (0, 1) and / = 1. In this case there holds 
u G ^^/■^^'"(r) for any r < 1/2 so that by Theorem 12.11 we expect a convergence of the mortar method 
close to /i^/^, the convergence of the conforming BEM, cf. Remark 1 2.11 This assumes that the mesh sizes 
h (of the sub-domain meshes) and k (of the meshes for the Lagrangian multiplier on the skeleton) are 
proportional, which will be the case in all our experiments. In fact, the elements of the mesh for the 
Lagrangian multiplier will always consist of two or three elements of the trace mesh. 

Since the exact solution u to (II. ip is unknown we approximate an upper bound for the semi-norm 
\u — Uh\ffi/2^q-y Here, we follow the strategy from [13]. Let us recall the procedure and discussion. 

By the ellipticity of V and (|3.8p there holds 

a{u -Uh,u- Uh) >\u- Uh\]ji/2(ry (^•^) 
On the other hand, using that u solves (jl.ip and (zVh solves (|2.ip . one finds 

a{u -Uh,u- Uh) = a{u, u) - 2a{u, Uh) + a{uh, Uh) = {Wu, u)r - 2a{u, Uh) + F{uh). 
By (13. 5p there holds 

a{u,Uh) = F{uh) - {[uh],\)^ 

such that, with the previous relation, 

a{u-Uh,u-Uh) = {Wu,u)r - F{uh) + 2{[uh]A)-f < {Wu,u)r - F{uh) + 2\\[uh]\\L^{^)\\\\\L^{^y). (4.2) 
Like in the proof of Theorem 13.21 one sees that ||A||/^2(^) is bounded. Therefore, by (14. ip we find that 

\u-Uh\]ji/2^r) - - F{uh)\ + \\[uh\\\L^{^y). 

The terms F{uh) and || [it/i] ||l2(^) are directly accessible and {Wu,u)y can be approximated by an ex- 
trapolated value that we denote by ||ti|jex (cf- [12]). Therefore, instead of the relative error 

Ik - '"fell//i/2(r)/ll^llHi/2(r)) 

we present results for the expression 

1 /2 

(l ll^llex ~ ^{'^h) \ + I|['"/i]||l2(7)) /II ^||ex 

(4.3) 

which is, up to a constant factor, an upper bound for \u — 'U/i|//i/2c2-)/|!^i||cx- 
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In the figures below we show different error curves, indicated by numbers (n) (n = 1, . . . , 4) as follows. 

(1) {\h\\L-F{uh)\+ WMWlH^))'^^ "mortar BEM" 

(2) {\\\u\\l,-F{uh)\y^^ "errorl" 

(3) IIKlll^/^) "error2" 

(4) a{u — Uh,u — UhY^"^ "conforming BEM" 

Here, Uh denotes a conforming boundary element solution. Additionally, all curves are normalized by 

1 1 ^ 1 1 ex • 

Therefore, to resume, an error curve (1) represents the upper bound (|4.3p for the (normalized) error 
1^ ~ ^/il_f/i/2(T) of the mortar BEM. Curves (2) and (3) are the two components of (1). Here, (3) controls 
the non-conformity of the mortar approximant u/j. Curve (4) represents the error of the conforming 
BEM. In this case it is equivalent to the error in energy norm ||u — 'U/i||j;^i/2(p^- 

All results are plotted on double logarithmic scales versus 1/h. For our numerical experiments we 
always use rectangular meshes and in this section, hi refers to the length of the longest edge on Ej, and 
h := maxj hi as before. 

Conforming sub-domain decomposition. 

Experiment 1 (conforming mesh, results in Figure [4?T]) . First let us consider a conforming decom- 
position of r into four sub-domains as indicated in Figure [2?TT a). Moreover, let us first test the case 
where the separate meshes on the sub-domains form globally conforming meshes (we take uniform meshes 
consisting of squares). The corresponding results are shown in Figure STJ Along with the curves (1), (2), 
(4) we plot the values of /i^/^. The numerical results indicate a convergence of the order 0(/i^/^), for the 
conforming as well as the mortar BEM. According to the discussion above this is the best one can expect. 
The curves (1) and (2), referring to our upper bound (|4.3p and the first term in (j4.3p . respectively, are 
almost identical. This means that the second term in (|4.3p . which in the next plots will be labeled by (3), 
is negligible in comparison. Indeed, in this symmetric case the jumps [uh\ disappear and the numerical 
results vanish at the order of single precision. Therefore, in this plot, we do not show the curve (3). 

We do not observe a logarithmical perturbation of the convergence in this range of number of un- 
knowns. This may be caused by the fact that we are not including the L^-parts in the error since our 
results are, up to constant factors, upper bounds only for the semi- norm \u — Uh\^i/2^q-y Also, we do not 
know whether our bounds including the logarithmic terms are sharp. 

Experiment 2 (non-conforming mesh, results in Figure [^?2]) . Now let us test globally non-conforming 
meshes. Again we use uniform meshes consisting of squares on each sub-domain. We mesh as in Fig- 
ure 12.1( a) starting with 2, 3, 4, and 5 "slides" on Fi, Fs, and respectively and increase the 
number of slides in each sub-domain by one in each step of our sequence of meshes. The corresponding 
results are shown in Figure l4?2l Again, a convergence of the expected order 0{h^/'^) is confirmed. Curve 
(3) indicates very fast convergence of the jumps || [^i/i] ||l2(^) — > 0. In the experiments below, however, 
we observe a slower convergence. In this particular sequence of meshes, where we increase the slides 
on the sub-domains by the same amount, the trace meshes from different sides on a particular interface 
edge approach each other in a certain sense. We conjecture that this specific situation ("approaching" 
conforming meshes) causes the fast convergence of the jumps. 

Experiment 3 (non-conforming mesh, results in Figure IT3]) . For the next experiment we start with 
a mesh of four squares on each sub-domain (the sub-domains are again as in Figure [2TTT a)). and increase 
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5 10 20 30 40 

1/h 

Figure 4.1: Conforming sub-domain decomposition with conforming meshes. 

the numbers of shdes on different sub-domains by different steps (increase by 2, 3, 4, 5 shdes on Fi, 
Fs, F4, respectively). In this case both error parts, curves (2) and (3), behave hke 0{h}/'^), confirming 
our a priori error estimate and thus the good performance of the mortar BEM. Let us note, however, 

1 /2 

that the part of the error expression (j4.3p is an overestimation. Indeed, our substitution ()4.3p 

for \u — 'W/i|/i-i/2('7-)/||u||ex is not precise. On the one hand we replaced the term 2 ||A||2,2(^) in ()4.2p by 1 
(and the generic constant in (j4.ip by 1). On the other hand the term ([u^], A)^ is of higher order than 
||[u/i]||^2(^). According to (I3.18[) and by standard approximation theory there holds for any r < 1/2 

L 

l(M,AKI < ||K]||i2(^) inf ||A-V'||l2(^) < r(j;|A|^.(^^))i/2||[^^j||^^^^^_ 

^ 1=1 

This shows that {[uh],X)^ is of higher order than ||[^i/i]||L2(^). Note that, by the proof of Theorem l3.2l and 
smce u E one has the regularity A G H/^i H'^ili) < 1/2- Therefore, by ()4.2p the term 

1 /2 

(l ||n||2^ - F{uh)\) /||n||ex (curve (2), "errorl") 

is asymptotically equal to 

a{u -Uh,U- M/t)^/^/||M||ex 

and this dominates the error. 
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Figure 4.2: Conforming sub-domain decomposition with non-conforming meshes, same refinement steps 
on sub-domains. 

Non-conforming sub-domain decomposition. 

Experiment 4 (non-conforming mesh, results in Figure . Finahy, we consider the fully non- 
conforming mortar method, i.e. with non-conforming sub-domain decomposition and non-conforming 
meshes. We decompose P into three sub-domains as in Figure and use the initial mesh given there on 
the left. Then slides on sub-domains are increased in each direction by 3, 2, 1 on Fi, F2, Fa, respectively, 
in each step. The second mesh is on the right in Figure Note that in each second step the cross-point 
(0, 0) between the sub-domains is a hanging node and our theory includes this case. The numerical results 
are shown in Figure 14.51 and again confirm the expected convergence of the mortar BEM. 

In this case, the meshes for the Lagrangian multiplier are coarsenings of the trace meshes from r2 on 
712 and 723, and of the trace mesh from Fi on 713. We always join two elements of the respective trace 
mesh to form an element of the Lagrangian multiplier mesh, except for an odd number of elements of 
the trace mesh when one set of three elements is joined. The corresponding numbers of unknowns for 
the steps are listed in Table [H 
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Figure 4.3: Conforming sub-domain decomposition with non-conforming meshes, different refinement 
steps on sub-domains. 



Figure 4.4: Conforming sub-domain decomposition with non-conforming meshes. 
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Figure 4.5: Non-conforming sub-domain decomposition with non-conforming meshes, different refinement 
steps on sub-domains. 



h = hs 


hi 


h2 


dim(X;,) 


dim(Mfc) 


0.5000 


0.1250 


0.1667 


27 


4 


0.3333 


0.0625 


0.0833 


80 


7 


0.2500 


0.0417 


0.0556 


161 


11 


0.2000 


0.0313 


0.0417 


270 


14 


0.1667 


0.0250 


0.0333 


407 


18 


0.1429 


0.0208 


0.0278 


572 


21 


0.1250 


0.0179 


0.0238 


765 


25 


0.1111 


0.0156 


0.0208 


986 


28 


0.1000 


0.0139 


0.0185 


1235 


32 


0.0909 


0.0125 


0.0167 


1512 


35 


0.0833 


0.0114 


0.0152 


1817 


39 


0.0769 


0.0104 


0.0139 


2150 


42 



Table 1: Dimensions and mesh sizes for experiment 4. 
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